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Abstract 

We propose a method for the determination of magnetic monopole currents in non- 
Abelian gauge theories which does not need a projection to Abelian degrees of freedom. 
With this definition we are able to determine the distribution of magnetic currents and elec- 
tric fields for the gluonic flux tube between a pair of static charges. Further we check the 



■ validity of the Gauss law and the dual London equation in a gauge invariant formulation. 



1 Introduction 



The success of the hypothesis of the existence of quarks and the failure of the detection of quarks 
as isolated objects is a big challenge for Quantum Chromodynamics. The formulation of QCD 
on a lattice provides the possibility to show that quarks are confined by a linearly rising potential 
p] the origin of which is the compression of the colourelectric flux between a pair of charges 
in a tube ||. The mechanism for the squeezing of the electric flux lines which originate in the 
colour charges is still a very interesting open question. One of the most promising conjectures 
is the hypothesis of 't Hooft @ and Mandelstam that the QCD vacuum behaves dually to a 
superconductor which expels magnetic flux lines by the Meissner effect. 

The validity of this dual superconductor picture can be tested much more easily in compact 
QED which is an Abelian simplification of the non- Abelian theory of QCD by using U (1) instead 
of SU(3) parallel transporters. In compact QED one can define magnetic monopole currents 
Q. Magnetic charges form a dual Cooper pair condensate in the confined phase and are partly 
expelled || by the dual Meissner effect from the region of the electric flux tube between a 
pair of electric charges. Curls of magnetic currents give rise to electric dipole moments. These 
electric dipoles align in the electric field of the charge pair. The electric field and the curl of 
the monopole currents are connected by the dual London equation Q. The aligned magnetic 
currents form a thin solenoid around the electric flux tube. 

In order to test the dual superconductor picture of confinement in non- Abelian gauge theories 
one has to invent a prescription for the identification of magnetic monopole currents. It was 't 
Hooft's idea || that the magnetic monopoles may be identified in a U (1) subgroup of the SU(N) 
gauge group. This leads to the question of an appropriate procedure for the determination of the 
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relevant degrees of freedom. Several partial gauge fixing procedures have been suggested which 
leave the largest Abelian subgroup unfixed. On the other hand a unique prescription is still 
missing. Moreover, in several articles problems appearing in the various gauge fixing procedures 
have been discussed ||. 

In order to avoid the gauge fixing problem we suggest in this article an alternative prescription 
|]i~0|| for the determination of magnetic monopole currents in non- Abelian gauge theories. This 
description is an extension of the method suggested in || for the identification of monopoles 
in f7(l) theory. We introduce our suggestion with a discussion of determination of the field 
strength and its connection to the Gauss law. Then we discuss the definition of gauge covariant 
monopole currents and finally we compare colourelectric fields and the curl of colourmagnetic 
monopole currents around a pair of colour charges and discuss the validity of the dual London 
equation. 



2 Field strength and the Gauss law 

In order to have a well settled definition of the field strength on the lattice we start our inves- 
tigation with the discussion of the Gauss law for a pair of static colour charges. In classical 
Chromodynamics the Gauss law can be derived by a variation of the full QCD action with 
respect to the forth component A^ a of the vector potential, with a being the colour index. To 
derive the Gauss law for Quantum Chromodynamics we use a version of Ehrenfest's theorem 
which is - for a gauge theory with gauge field U - given by 

with O being an arbitrary operator. In the path integral formalism (p]) is a consequence of the 
translation invariance of the integration measure of the path integral. On a four-dimensional 
Euclidean hypercubic lattice A of spacing a with periodic boundary conditions in space and time 
direction the gauge fields are defined on links (x, fj) by parallel transporters 

U^x) = e ia JW F " F a = ^, a%(x) = agAftx), x G A, M = 1, 2, 3, 4 (2) 

where A a denotes the Gell-Mann matrices and g the QCD coupling constant. An action conve- 
nient for our purpose is the standard Wilson action 

SP] = P £ (l-^ReTr[V(x)) . P = "| ( 3 ) 

x^<v \ 6 9 



with U fJjU (x) being the product of link variables around an elementary plaquette in /i^-direction 
at lattice site x 

U^(x) = U^{x)U v {x + fi)Ul{x + i))Ul{x). (4) 

Since we consider pure SU(3) gauge theory with external static colour charges, the operator O 
in (H) is given by the Polyakov loop 

L(x+) = J] U 4 (x+,t) (5) 
t=i 
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describing the time evolution of a static charge at lattice site x+. To derive the Gauss law at 
an arbitrary point x the Polyakov loop (||) has to be parallel transported to x by a Schwinger 
line U(x,x+) and reads 

L x {x + ) = U(x,x + )L(x+)U(x+,x), (6) 

with U(x, x + ) being the product of link variables along a path connecting the lattice sites x and 
x+. For the variation in (|l|) we have to write 

g - 9 X (7\ 

5a%{x') dai(x) M '" [ } 

Summation in colour space guarantees a gauge invariant formulation of the Gauss law and we 
receive for x/x- for (jl]) the following expression 

^]^^«>^x + ))TVL*(s_),r^ = 

= Jv{U]- j ^- ) Tv(F a L*(x + ))TTL'(x-)e- s < u 'l. (8) 

The case x = X- is treated in section 4. The approximation 

leads to 

iCf f) J V[U] TrL(x+) TrL*(x_) e" 5 ^ 1 5 X , X+ (10) 
for the right-hand side of expression (|8|), with 

Cf f) := F a F a = i, (11) 

being the eigenvalue of the quadratic Casimir operator in the fundamental representation. To 
find a simplification for the left-hand side of (|8|) we keep in mind that a link variable U±{x) 
contributes to six space-time plaquettes U^. This we have to consider carrying out the derivative 
of the action S. Using 

F a Tr(F a U^(x)) = itV(x) (tJ) (12) 
we obtain the following result for the Gauss law (|8]) 

( Tr | \ga 2 D'wE(x)\ X L x {x + ) \ TrL*(f_) ) 4 

(TrL(x + )TrL*(x_)) 3 '+ r V ; 

where we introduced the definition 

(DivE(x)) X = X; (^(x) - Ef(x - ij) (14) 
i=l 

for the covariant divergence of the electric field strength. The electric field strength itself turns 
out to be the following hermitean traceless operator 

ga 2 EUx) = j. (U u (x) - U^x)) (15) 
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Figure 1: The worldlines of the 12 contributions in TV < ga Div£'(x+) L(x+) \- 



and is directly connected with the action (||) used in our derivation of the Gauss law. If we 
want to know the field strength Ei(x) at the lattice site x + i, we have to carry out a parallel 
transport to the corresponding site which reads 



E* + \x) = Uj(x)Ef(x)U l (x) 



(16) 



and is nontrivial for a non Abelian gauge theory. Of course, the field strength ([D]) is an element 
of the s«(3)-algebra and therefore as an 8-vector in the algebra a gauge dependent quantity. 
But the correlation ( |l3|) shows us how to calculate it as a gauge invariant quantity in lattice 

ga 2 Y)i\E(x + ) + in ( |l3| ) can be represented by the sum of 

Tie formulation of the Gauss law (|i"3| ) appears to be asymmetric with 
respect to the positions of quark and antiquark. According to our interpretation the Polyakov 
loop at x+ defines a direction in the eight-dimensional group space. The average projection of 
T)ivE(x) to this direction is non zero and given by C 2 « for x = x + only. In other words the 
Polyakov loop is the source of an electric field which points in average in the direction given by 



simulations. The term Tr j 
worldlines shown in FigfiJ. 



the Polyakov loop itself. The antiquark in (13) is only a "spectator" which is the necessary sink 
of the flux lines originating in x+. 

The right-hand side of the Gauss law (|l3|) corresponds to the square of the charge density of 
the considered colour sources. The charge density turns out to be imaginary as it should be in 
Euclidean space|] |l^] and is proportional to the eigenvalue of the quadratic Casimir operator 



= | of a colour triplet which is a direct consequence of carrying out the trace in colour 
space in (|8|). 

The Gauss law can also be derived for a static colour octet. We consider a 3 
lattice site x with generator 

Fj}®~ 3 = F a ® 1 + 1 (2) F a 



and the Polyakov loops for this composed system 

L 3 ® 3 (x) = L{x)®L*{x) 



3 -system at 
(17) 

(18) 



f We use the following relation between Minkowski and Euclidean observables: All time components of 4- vectors 
transform like x 4 = ix . In addition we use the convention (E, B, p m , J m ) Euclidean = (iE, B, p m , —iJ m ) Minkowski- 
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and obtain in analogy to 
V[U] Tr 



a da%{x) K ' 



-S[U] 



V[U]Tr 



F 



3®3 







L 



3®3 1 



-S[U] 



(19) 



a da%{x 

Since the uncharged singlet does not contribute, this leads to the Gauss law for the colour octet 
( 2i 9 {Tr [(ga 2 DivE(x)) x L(x)] TrL*(x)} 



TrL(x)TiL*(x) - 1 } 



i3g 2 . 



(20) 



As expected the colour charge on the right-hand side of (|20| ) is proportional to the eigenvalue of 
the quadratic Casimir operator Cf a \ = 3 in the adjoint representation. 



3 Definition of magnetic monopoles 

To identify magnetic monopoles in Abelian gauge theories DeGrand and Toussaint || suggested 
to count the number of Dirac strings emanating from a three-dimensional cube. The Dirac 
strings carry a quantized magnetic flux of — . The net flux for a closed surface defines the 
number of magnetic monopoles within a cube as a gauge invariant quantity, 't Hooft suggested 
to identify magnetic monopoles in non- Abelian gauge theories by a gauge fixing procedure where 
an appropriate operator in the adjoint representation of the gauge group is diagonalised. The 
monopoles are then defined as Dirac monopoles of the remaining [/(l)-symmetry. The problem 
of this formalism is that the monopoles depend on the procedure used to define the Abelian 
projection. 

In H an alternative definition of monopoles in Abelian gauge theories was suggested which 
uses the dual Maxwell equations. The monopole currents J m n in this definition read on the 
lattice 

ZfiupcrO 2 51 = —Jm,fi, ^1234 = +1 (21) 

□ i£cube(i/p<r) 

with i ? n i being the field strength according to the definition 

a 2 eF^ = sin^ (22) 

where a is the lattice spacing, e the electric coupling and 9^ the angle of an elementary plaquette 
in the /iz^-plane. In the following we describe a generalization of ( |2~l| ) to non-Abelian gauge 
theories as QCD. According to ( |i~5|) we define 

ga 2 F* v {x) = i(lV(x)-[/t,(x)) (H) , = 1,2,3,4 (23) 

i.e. 

Ef (x) = F%(x), Bf(x) = ^e ijk Ff k {x). i,j,k = 1,2,3 (24) 
With the covariant derivative of the field strength 

D v F%{x) = F*{x + 0)- F*(x) = U u {x)F*+\x + u)Ut{x) - F%{x). (25) 
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Figure 2: The types of paths contributing to a three-dimensional spatial monopole cube, defining 
together with the reversed paths according to (23) the monopole current in time direction. 



the dual Maxwell equations in generalisation to (21) read 

1 



2 



e^D^x) = -J^(x). (26) 



They define the colourmagnetic monopole current J m . M . On the lattice each component of J m ^ 
corresponds to a three-dimensional cube built of six plaquettes which measure the flux out of 
the cube, see Fig.|2[ The covariant derivative ( |25|) contains a parallel transport and therefore 
guarantees that all contributions to the total flux out of the cube are determined at the same 
lattice site x. 

According to the Bianchi identities the magnetic currents ( |26| ) should identically vanish. The 
lattice results reported below show that for finite lattice spacing this is obviously not the case. 
In the continuum the Bianchi identities are derived from the definition of the field strength 
via potentials. Also on the lattice they should be regained for "infinitely small" cubes. There 
may be two reasons why the magnetic currents do not vanish identically in lattice calculations: 
Bianchi identities may be violated by Abelian monopoles and non-Abelian monopoles. Abelian 
monopoles are attached to Dirac strings, i.e. to 2-7r-rotations of Abelian potentials. Such Abelian 
monopoles could appear in an arbitrary [/(l)-subgroup of the full gauge group. They are point- 
like and could be regarded as lattice artifacts. The non-Abelian monopoles could be associated 
with the non-Abelian nature of the gauge fields. To make this plausible we cover a cube by three 
pairs of plaquettes which are surrounded by paths as shown in Fig.||. For the parallel transport 
along path Cj we write 

Ui = e m , Ti € su{3), i = 1,2,3. (27) 

It is obvious that 

e m e iF 2e iF 3 = 1 (28) 
Nevertheless, according to the Baker-Campbell-Hausdorff formula we have in general 

Tx + T 2 + + 0, (29) 
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Figure 3: Three pairs of plaquettes covering a monopole cube and measuring the flux out of the 
cube. 

what leads to a non-vanishing contribution to the magnetic current which should vanish for 
infinitely small cubes. Therefore, such non-Abelian monopoles exist in the continuum limit only 
if they are extended objects. 

It is an interesting question if the current defined in (^6[) is conserved. The flux through a given 




Figure 4: (a) One plaquette contributes to four monopole cubes or - seen from another point of 
view - the four monopole cubes contributing to a given component of the curl have one plaquette 
in common, (b) The eight monopole "cubes" contributing to the four-dimensional divergence of 
the monopole current on the lattice. Since every plaquette contributes twice, once with positive, 
once with negative sign, the current is conserved. 

plaquette contributes to four adherent monopole cubes, see Fig.||a, which form the links of a 
dual plaquette on a dual lattice. The corresponding contribution to the magnetic current runs 
cyclically around this dual plaquette and is therefore conserved. We can discuss this problem 
also in another way. On the lattice there are eight magnetic current components contributing 
to the continuity equation. The corresponding eight cubes are depicted in Fig.|]b, where the 
time evolution can be identified with a magnification of the small cube to the large cube. It 
can be easily seen that every plaquette contributes to two cubes with positive and negative sign 
respectively. If one guarantees that both contributions of a given plaquette are transported on 
the same path to the given reference point the sum of the 48 contributions is zero. 
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Since our aim is to verify the validity of the dual London equation, we need to define not 
only the colourmagnetic monopole currents J mjV (26) but also the curl of the monopole current. 
We use the following discretisation of the curl 

(curlJ™) 2 (x) = e ijk (DjJ mtk ) x = e ijk (j^ k {x) - J^ k (x - j)j i,j,k = 1,2,3 (30) 

with Dj being the covariant derivative as defined in (^), but we simultaneously take into account 
that the field strength contributions to the current on the right-hand side of (30) have to be 
transported along the shortest path to the reference site x of the curl. For the further discussion 
of the London equation we want to emphasize that on the lattice the four monopole cubes 
contributing according to ( |30| ) to a given component of the curl have one plaquette in common, 
as shown in Fig.||a. This leads to a correspondence of the operators (curlJ m )f (a;) < — > Ef(x) 
which refer both to the same lattice site x. 

We have now defined all necessary ingredients for a numerical determination of the electric 
field distribution and the distribution of the curl of the magnetic current and to test in the 
following section the validity of the dual London equation for the gluonic vacuum with external 
static colour charges. 



4 Correlation functions and numerical results 

In this section we would like to present our numerical results. All simulations were done on a 
Euclidean 8 3 x 2-lattice in the confinement phase at (3 = 4.9. The gauge field configurations 
were generated with a Hybrid-Monte-Carlo-algorithm. 

First let us begin with the "measurements" concerning the verification of the Gauss law for a 
colour triplet respectively a colour octet. The left diagram of Fig.|| shows the results of evaluating 
the correlation function ( |l3| ) of 40.000 gauge field configurations for different distances R[a] of 



\(ga 2 T)\vE{x)) triplet + 

2.5 - 



2 - 

+--+--+-+-+-+--+-+-+++-+-+- 
1.5 - 



1 Li 

1 2 3 4 10 20 30 40 

R[a] Number of trajectories in 1000 




Figure 5: The Gauss law for a colour triplet as a function of the distance R[a] between the 
charges (left figure, see Eq. (13)) and for a colour octet as a function of the number of evaluated 
gauge field configurations (right figure, see Eq. (20)). In both cases the theoretical prediction 
is plotted as a dotted line and agrees very well with the numerical values. 

the colour charges. The theoretical prediction is plotted as a dotted line and agrees very well 
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Figure 6: A Schwinger line U(x+, x) connects the field strength Ei(x) = ^ (u^x) — U^x) 
and the Polyakov loop L(x+) to guarantee gauge invariance. 



(ti) 



with the numerical results of the left-hand side of The right diagram of Fig.|5] shows the 

numerical value of the correlation function ( p0|) as a function of the number of evaluated gauge 
field configurations. The results agree with the theoretical value. 

Now let's turn to the correlation functions measuring the distributions of the electric field and 
the curl of the the monopole current. The derivation of the Gauss law on the lattice demonstrates 
how to define the electric field strength as an 8- vector in colour space (|l5|). To determine it as 
a gauge invariant quantity we have to evaluate - according to (|l3|) - the correlation function 



QQ 



(Tr (u(x + ,x)Ef(x)W{x+,x)L(x + )^ TrL*(x_)) 



(TrL(x+)TrL*(£_)) 



(31) 



with U(x+,x) being the Schwinger line connecting the field strength Ef{x) with the Polyakov 
loop L{x + ) (Fig.|6]). This correlation is very similar to the one proposed in |TT| for the electric 
field. To measure the distribution of the curl of the monopole current we generalize (^) and 
replace the electric field strength by the curl defined in pC ) . Because of local gauge invariance of 
QCD in colour space we have to refer both operators in (3l|), Polyakov loop and field strength, 
to the same lattice site. This is guaranteed by the Schwinger line in (pl|). But there is no unique 
choice of Schwinger line connecting the two operators. Since we try to verify the validity of 
the dual London equation what means comparing two quantities with each other, this problem 
has no consequences to our intention, if the electric field strength and the curl of the monopole 
current are transported on the same path and therefore measured in the same local coordinate 
system. There is another property of the correlation ( |3l| ) which should be mentioned. Charge 
and anticharge do not play a symmetric role as one would naively expect, because the electric 
field can only be connected with either the charge or the anticharge. Therefore we should not 
expect a result which is symmetric with respect to the two charges. 

Now we consider the numerical results which we obtain from evaluating the correlation 
function (|3l| ) for the electric field and the curl of the monopole current. We used 50.000 gauge 
field configurations and varied the distance between the charges from la to 4a. In Figj7| we 
show a transverse profile of the electric field strength (left column of Fig.0) and the curl of the 
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Figure 7: Transverse profile of the electric field strength (left column) and the curl of the 
monopole current (right column) in the plane perpendicular to the QQ-axis for charge distances 
d = la, d = 2a and d = 3a. The electric field strength decreases quickly with increasing radial 
distance from the QQ-axis. This indicates that the flux tube is mainly concentrated on the 
link(s) connecting the charges. The absolute value of the curl shows the same behaviour as 
the electric field strength, only on the link(s) connecting the charges they have a different sign 
because of the fluxoid. The errorbars shown indicate the jack knife error, for d = la and d = 2a 
they have been omitted, because they are smaller than the plotting symbols. 
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Figure 8: Longitudinal profile of the electric field strength (left column) and the curl of the 
monopole current (right column) on the QQ-axis for charge distances d = la, d = 2a and 
d = 3a. For d = la the charges are at zeta = ±0.5a, for d = 2a at zeta = ±la and for d = 3a at 
zeta = ±1.5a. For all three charge distances the electric flux as well as the curl of the monopole 
current is mainly concentrated on the link(s) connecting the two charges. As described in the 
text the full (dotted) line corresponds to a symmetrized (non symmetrized) "measurement". 
The error bars shown indicate the jack knife error. 
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monopole current (right column of Fig.]?]) in the plane perpendicular to the axis connecting the 
charges for charge distances d = la, d = 2a and d = 3a. The electric field strength decreases 
quickly with increasing radial distance from the QQ-axis what characterizes a thin flux tube 
between the charges. The absolute value of the curl of the monopole current shows the same 
behaviour as the electric field, only on the QQ-axis the curl is reversed due to the fluxoid 0. 
In Fig|| we show the electric field strength (left column of Fig.||) and the curl of the monopole 
current (right column of FigJ^) on the QQ-axis for charge distances d = la, d = 2a and d = 3a. 
The full line represents a symmetrized "measurement", where the electric field strength (the 
curl of the monopole current) is connected to the close charge. The dotted line depicts the 
result for Schwinger lines connecting the field strength (the curl of the monopole current) with 
the far charge. As expected the sign of the electric field strength is changing at the position of 
charges. In the case d = la the flux is confined almost entirely within the link connecting the 
two charges. Only 2% of the total flux reach the anticharge on the long way through the lattice 
boundary. For a charge distance of d = 3a the electric flux as well as the curl of the monopole 
current decreases in the middle between the charges what can be interpreted as fluctuations of 
the fluxoid. We want to emphasize that the curl of the monopole current has the same sign 
along the whole axis; according to the definition (|26| ) it forms a solenoid of left handed currents 
after going back to Minkowski space. 

As we discussed in connection with the Gauss law a triplet or octet Polyakov line is the source 
of a coulor electric flux in a certain direction in the eight-dimensional su(3)-algebra space. The 
electric flux is continuously "rotated" on its path to the anticharge due to the non Abelian 
nature of the group SU(3). A demonstration of this effect is the lowering of the dotted line 
compared to the full line in Fig.|8] showing the profile of the electric field and the curl of the 
monopole current along the QQ-axis. Another evidence for this effect can be obtained from the 
evaluation of (|8|) for x = x_ resulting in 

(Tr{[oa 2 Div£(x_)P L{x + )) TtL*(x _) ) = -i± g 2 (Tr [l x ~ (x + ) (tl) L\x. 

= -i X -g 2 ( Tr [L*-(x+)tf(x-j\ - 1 -TtL(x + )TtL*(x^) ) 

= -i\g 2 {3[L x -{x+)®L*{xJ)] singlet - ~Tr£(x + )TrL*(x_) ) 

= -i\ 9 2 { [L*-(x + )®L*(x-)] singlet - ±Tr[L*-(x + )®L*(x-)] octet ). (32) 

The last expression turns out to be proportional to the difference of Wilson loops in the singlet 
and octet representation of a quark-antiquark system at distance d =| x + — x_ |. In Fig.^ we 
show the numerical results for the absolute value of the left hand side of (^) divided by 

<TrL(x + )TrL*(x_) ) = ( [L*-{x+) L* {x„)] singlet + ^Tr [L*-(x+) ® L*(x-)] octet ). (33) 

The divergence produced by the anticharge, the sink of the electric field, is an s-u(3)-vector which 
with increasing distance d =| x+ — x_ | loses its correlation with the direction defined by the 
Polyakov line matrix at x + , the source of the electric field. 



The connection between DivE(x^) x+ and the QQ potential in (32) together with the correlation 
of two Polyakov lines in (|33|) gives the nice opportunity to extract the potential of a QQ-pair in 
the singlet and octet channel 
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Figure 9: The divergence of the electric field at the positions of a quark and an antiquark in 
dependence on the distance between the charges. In both cases the gauge dependent divergence 
is connected to the Polyakov line of the positive charge by a Schwinger line in order to define a 
gauge invariant quantity. 



V s inglet(\ X+ - X- I) 
V oc tet(\ X+-X- I) 



1 



— ln([L x -(x+)®L*(x-)] singlet ), 

— ln(lTr[L x -(x + )^L*(x^)] octet ), 



N t a 



(34) 



with Nt being the number of lattice sites in time direction. Because of the periodic boundary 
conditions on the 8 3 x 2-lattice the QQ-pair of (|33| ) may be connected by two different strings 
of length ri =| x + — x_ [ and T2 = N x a — r±, where N x denotes the number of lattice sites 
in space direction. For charge distances smaller than N x a/2 we can approximately neglect the 
contribution of the longer distance. For the distance N x a/2 we have to consider only half of the 
value of the left-hand side of ( |33|) for the determination of the potentials in (|34|). The numerical 
results of are shown in in Fig. JO]. 

Since our aim is to test the dual superconductor picture of confinement, we turn now to the 
dual version of the London equation. Taking into account fluxoid contributions it reads in 
Minkowski spaceQ 



{E)qq = A 2 (curlJ m )QQ + fluxoid contributions 



(35) 



where A is the London penetration depth. In Fig.|ll| and in Fig. 12 we show the ratio of the 
electric field strength and the curl of the monopole current on the links of the three-dimensional 
lattice for charge distances d = la and d = 2a, where only results from symmetrized calculations 
are considered. It is clearly seen that the ratio 



A 



\^>QQ/(( curlJ « 



IQQ 



(36) 



^According to our convention for the relation between Minkowski and Euclidean observables, in Euclidean 
space the dual London equation reads (E)qq = —A 2 (curlj m )gg. 
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Figure 10: QQ-p&ir potential in the singlet and octet channel as a function of the distance 
between the charges. 



rises with the distance d. (Be aware of the different scale in Fig.|TI| and in FigjT^.) The same 
increase for A 2 has been found in compact QED Q and can be interpreted as increasing effect of 
string fluctuations. For the case d = la, off axis A 2 lies in the range A 2 = 0.08 ±0.02. This leads 
to a fluxoid on the link between the charges $ = iga 2 ({Ei) qq — A 2 ((curlJ m )j) qq) = 1.604±0.077 
which agrees quite well with the theoretical value of <& = | g 2 = 1.633. For a charge distance 
d = 2a we obtain the value A 2 = 0.14 ± 0.03. The value of the corresponding fluxoid is then 
given by $ = 1.504 ± 0.106. A fit in an effective model allowing string fluctuations has shown 
that in compact QED |j] the London penetration depth A is even smaller than A for distance 
d = la. The above results show that the ratio (E^qq/ '((curlJ m )j)QQ behaves analogously to 
compact QED. According to our opinion this demonstrates the similarity of the confinement 
mechanism in QCD and compact QED. 



5 Conclusion 

In this article we developed a formalism how to determine the electric field distribution and the 
distribution of the curl of the monopole current of the gluonic vacuum with an external QQ-p&ii 
without using the technique of Abelian projection. We started with a derivation of the Gauss 
law for a colour triplet to find a well settled definition for the field strength on the lattice. By 
using the dual Maxwell equations we defined colourmagnetic monopoles and further the curl 
of the monopole current. We introduced Schwinger lines to determine the colourelectric field 
strength and the curl of the monopole current as gauge invariant quantities in lattice simulations 
and found that the ratio (E^qq/ ( (cur \Jm)i)QQ behaves as predicted by the dual superconductor 
picture of confinement. 

In future investigations we would like to measure the density of colourmagnetic monopoles 
defined in (|26|). A candidate for a gauge invariant operator is the length of the monopole current 
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Figure 11: Ratio of the electric field strength and the curl of the monopole current on the links 
of the three-dimensional lattice for charge distance d = la. 
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Figure 12: Ratio of the electric field strength and the curl of the monopole current on the links 
of the three-dimensional lattice for charge distance d = 2a. 



15 



in su(3)-space | J M | = yJ2a=i Jjl^jl- First calculations show that the value of (| |) decreases 
with increasing inverse coupling (5 crossing the phase transition. But the quantity | | has the 
disadvantage that there are only positive contributions. Quantum fluctuations which are not of 
topological origin do not cancel. We want to point out that for the correlation function (j3l|) 
the just mentioned problem does not occur: The Polyakov line fixes a direction in colour space, 
fluctuations of the operator of the electric field strength or the curl of the monopole current 
contribute equally with positive and negative sign and therefore cancel, whereas contributions 
of topological origin survive the averaging process. 
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